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1. Introduction 


Wireless networked systems arise in various communication contexts, and are becoming 
a bigger and integral part of our everyday life. In today practical networked 
systems, information delivery is accomplished through routing: network nodes simply 
store-and-forward data, and processing is accomplished only at the end nodes. Network 
Coding (NC) is a recent field in electrical engineering and computer science that breaks 
with this assumption: instead of simply forwarding data, intermediate network nodes may 
recombine several input packets into one or several output packets (Ahlswede et al., 2000). NC 
offers the promise of improved performance over conventional network routing techniques. 
In particular, NC principles can significantly impact the next-generation of wireless ad hoc, 
sensor, and cellular networks, in terms of both energy efficiency and throughput (Ho et al., 
2003). 

However, besides the many potential advantages and applications of NC over classical 
routing, the NC principle is not without its drawbacks. A fundamental problem that NC 
needs to face over lossy (e.g., wireless) networks is the so-called error propagation problem: 
corrupted packets injected by some intermediate nodes might propagate through the network 
until the destination, and might render impossible to decode the original information (Cai & 
Yeung, 2002). As a matter of fact, the application of NC to a wireless context needs to take into 
account that the wireless medium is highly unpredictable and inhospitable for adopting the 
existing NC algorithms, which have mostly been designed by assuming wired (i.e., error-free) 
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networks as the blueprint. Furthermore, in contrast to routing, this problem is crucial in 
NC due to the algebraic operations performed by the internal nodes of the network: the 
mixing of packets within the network makes every packet flowing through it statistically 
dependent on other packets, so that even a single erroneous packet might affect the correct 
detection of all the other packets (Koetter & Kschischang, 2008). On the contrary, the same 
error in networks using just routing would affect only a single source-to-destination path. An 
up-to-date overview of recent results and open problems related to the application of NC toa 
wireless context can be found in (Di Renzo et al., 2010a), (Di Renzo et al., 2010b). 

NC finds successful application in wireless cooperative networks (Pabst, 2004), (Scaglione et 
al., 2006), since it offers an efficient way for overcoming their limitations in terms of achievable 
throughput (Katti et al., 2008b). More specifically, in wireless cooperative networks multiple 
radios are deployed in a neighborhood. The radios connect to each other via wireless links 
to form a multi-hop wireless network, with a few nodes acting as gateways that connect the 
wireless network to, e.g., the Internet. Packets traverse multiple wireless links before reaching 
the gateway and finally the wired network (and, thus, the destination). Multi-hop networks 
extend the coverage area without expensive wiring, thus offering cheap and moderately fast 
connectivity, which is often sufficient for accessing the Internet assuming normal browsing 
habits (Pabst, 2004). However, the price to be paid for this improved robustness in the 
transmission of data is the reduction of the achievable data rate, which is due to two main 
reasons: i) the transmission of redundant information to achieve the benefits of spatial 
diversity, and ii) the practical need to adopt the half-duplex constraint, which precludes the 
nodes to transmit and receive data simultaneously (Wang et al., 2009). In this context, NC 
offers an intelligent solution to boost the channel capacity of multi-hop wireless networks by 
combining information from different packets at the packet (Katti et al., 2008a), (Chachulski 
et al., 2007), symbol (Katti et al., 2008c), or signal (Katti et al., 2007) level: data mixing allows 
the system to offset the throughput limitations set by the half-duplex constraint and to reduce 
the amount of redundant data to be transmitted. For example, the inherent capability offered 
by NC to recover the throughput of cooperative networking has been clearly assessed for the 
so-called two-way relay channel (Zhang et al., 2006). 

Moving from the considerations above, it is evident that the design of efficient and robust 
protocols and algorithms to exploit the properties of NC in a cooperative networking scenario 
will play an important role for the next generation of wireless networks that require high 
data transmission rates. In particular, to take full advantage of the benefits of NC in a 
cooperative scenario, the network codes have to be properly designed in order to: i) maximize 
the probability of correct decoding at the destination nodes in order to minimize the effects of 
the error propagation problem, ii) reduce the energy consumptions of the overall cooperative 
network with the aim to prolong its operational life, and iii) keep the complexity of the 
relay nodes performing NC at a low level. Originally, the design of network codes has 
mainly been concerned with methods to achieve the maximum information flow (Ahlswede 
et al., 2000), (Li et al., 2003), (Koetter & Medard, 2003), (Ho et al., 2006). However, in 
the recent period considerable effort has been devoted to the design of efficient network 
codes to attain the maximum diversity gain (Xiao & Skoglund, 2009a), (Xiao & Skoglund, 
2009b), Rebelatto et al. (2010a), (Rebelatto et al., 2010b), (Topakkaya & Wang, 2010), which 
is known to determine the Bit Error Probability (BEP) for high Signal-to-Noise-Ratios (SNRs) 
(Wang & Giannakis, 2003). More specifically, as far as a multi-source multi-relay cooperative 
scenario is concerned, in (Xiao & Skoglund, 2009a) it has been shown that binary NC is 
sub-optimal for achieving full-diversity, and in (Rebelatto et al., 2010b) it has been pointed 
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out that max-diversity-achieving network codes can be obtained by resorting to the theory 
of non-binary linear block codes. For example, as far as the canonical two-source two-relay 
cooperative network is concerned, the methods proposed in (Xiao & Skoglund, 2009a), 
(Xiao & Skoglund, 2009b), (Rebelatto et al., 2010a), (Rebelatto et al., 2010b) can achieve 
full-diversity, when, instead, XOR-based binary NC (Katti et al., 2008a) cannot. The solution 
proposed by all these papers to overcome the limitation in the achievable diversity is based 
on using network codes in a non-binary Galois field. However, the price to be paid for 
this performance improvement is the additional complexity required at the relay nodes, 
which must network-code the received packets by using non-binary arithmetic. Also, longer 
decoding delays are, in general, required to design full-diversity-achieving network codes 
(Rebelatto et al., 2010b). Furthermore, the solutions available in (Xiao & Skoglund, 2009a), 
(Xiao & Skoglund, 2009b), Rebelatto et al. (2010a), (Rebelatto et al., 2010b) aim at guaranteeing 
the same diversity gain for all the active sources of the network, which, in general, leads to 
an inflexible network code design as multiple sources might have different Quality-of-Service 
(QoS) requirements, and, so, might need different diversity gains. In conclusion, the solutions 
available so far seem to be still inflexible to accommodate the needs of the multi-source 
multi-relay scenario, as well as to keep the computational complexity of the relay nodes at 
a low level. 

Motivated by these design challenges, the aim of this book chapter is to propose a new and 
flexible method to design network codes for cooperative wireless networks with the objectives 
of: i) improving the diversity gain of conventional relay-only (Scaglione et al., 2006) and 
XOR-based binary NC (Katti et al., 2008a), ii) keeping at a low level the complexity of the 
relays, and iii) having the flexibility of assigning to each source a different diversity gain 
according to the desired QoS. In this book chapter, we show that these design goals can be 
simultaneously achieved by exploiting the theory of Unequal Error Protection (UEP) linear 
block codes for the flexible and robust design of network codes for multi-source multi-relay 
networks (Masnick & Wolf, 1967), (Boyarinov & Katsman, 1981). In particular, by focusing 
on the canonical two-source two-relay network scenario we prove that, by adopting a simple 
(4,2,2) UEP code (Van Gils, 1983, Table I) as a network code, at least one source can achieve 
a better diversity gain than conventional relay-only or XOR-only NC protocols without the 
need to either use non-binary operations or require additional time-slots. The adoption of 
UEP coding theory for wireless relay networks (Nguyen et al., 2010) and for random NC 
with application to multimedia content distribution (Thomos & Frossard, 2004) is receiving a 
growing interest. However, to the best of the authors knowledge, none of the available works 
have exploited UEP coding theory for the flexible design of distributed network codes for 
diversity purposes. 

The reminder of this book chapter is organized as follows. In Section 2 and Section 3, system 
model and network code design are introduced, respectively. In Section 4, a low-complexity 
detector based on the Maximum-Likelihood (ML) principle for demodulation and network 
decoding at the destination node is derived. In Section 5, an analytical framework to compute 
the Average BEP (ABEP) and the diversity gain of various network codes is proposed. In 
Section 6, numerical results are shown to substantiate our claims and compare the UEP-based 
network code design with conventional relay and NC methods. Finally, Section 7 concludes 
the book chapter. 
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Fig. 1. Two-source two-relay network topology. Different line-styles denote transmission 
over orthogonal channels (e.g., time-slots (Scaglione et al., 2006)) to avoid mutual 
interference: S4 transmits in time-slot 1 (solid lines), S5 in time-slot 2 (dashed lines), R4 in 
time-slot 3 (dotted lines), and R» in time-slot 4 (dashed-dotted lines). 


2. System model 


Let us consider the canonical two-source two-relay cooperative network shown in Figure 1. 
Generalization to multi-source multi-relay networks is possible, but it is not considered in 
this book chapter due to space constraints. The working principle of the network in Figure 1 
is as follows. In time-slots t = 1, 2, source node S; broadcasts a modulated symbol, xs,, with 
average energy Em. For analytical tractability, Binary Phase Shift Keying (BPSK) modulation 
is considered, i.e.: 


xs, = VEm(1 — 2bs,) (1) 
where bs, = {0,1} is the bit emitted by S;. 


Accordingly, the signals received at relays R1, R?, and destination D are: 


Vs,R, = ASR Xs, + NSR 
Vs,R; = ASR XS, + fts,R; (2) 
Vs,D = ls,pXs, + is, 


where /1xy is the fading coefficient from node X to node Y, which is a circular symmetric 
complex Gaussian random variable with zero mean and variance 0% per dimension (i.e., a 
Rayleigh fading channel model is considered). For analytical tractability, independent and 
identically distributed (i.i.d.) fading over all the wireless links is considered, i.e., o5 = gy 
for any X and Y. Furthermore, n xy denotes the complex Additive White Gaussian Noise 
(AWGN) at the input of node Y and related to the transmission from node X to node Y. The 
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AWGNs in different time-slots are independent and identically distributed with zero mean 
and variance No /2 per dimension. 

Upon reception of ys,n,, Vs,R,, Vs,R,, aNd Ys,R,, the relays R and Rz attempt to decode 
the symbols transmitted by Sı and Sz in a similar fashion as in a Decode-and-Forward 
(DF) cooperative protocol (Scaglione et al., 2006). Unlike other solutions available in the 
literature for network code design for cooperative networks (Xiao & Skoglund, 2009a) (Xiao & 
Skoglund, 2009b), (Rebelatto et al., 20102), (Rebelatto et al., 2010b), we do not rely on powerful 
(i.e., Shannon-like) channel codes at the physical layer, which allow each relay to detect correct 
and wrong packets, and enable them to forward only the former ones. We consider a very 
simple implementation in which the relays demodulate-network-code-and-forward (D-NC-F) 
each received symbol without checking whether the symbol is correct or wrong. The main 
aim of this assumption is to keep the complexity of the relays at a very low level (Koetter & 
Kschischang, 2008), and to understand the robustness, in terms of coding and diversity gain, 
of the error propagation problem (Di Renzo et al., 2010a), (Di Renzo et al., 2010b) on various 
network codes. 

In particular, we consider ML-optimum demodulation that exploits Channel State 
Information (CSI) about the source-to-relay channels at each relay node, as follows (t — 1, 2): 


bs,R, = argmin( |ys,n, =A Emhs,r, (1 = 2bs,)|7} 


bs,e{0,1} (3) 
Îs,R, = argmin{|¥s,R, ~ yV Emhs,r,(1 EE 2bs,) |? 
bs,e{0,1} 


where ^ denotes the estimated symbol, and ~ denotes the trial symbol used in the 
hypothesis-detection problem. 
After demodulating bs, for t = 1,2 and q = 1,2, each relay, Rg, executes the following 
operations: i) it performs NC on these symbols, ii) it re-modulates the network-coded symbol, 
and iii) it transmits the modulated symbol to the destination during the third (if q = 1) and 
fourth (if q = 2) time-slot. By denoting with fr, (-,-) the NC operation performed by relay Ry, 
and with b R, the network-coded symbol, i.e., b R, = f. R, (bs, Ry bs, R)» the signal received at the 
destination, D, is: 

VR,D = hn,DXR, + RD (4) 


where xg, — V Em (1— 2br, ). 


After four time-slots, the destination has four received signals, i.e., Vs, p, Vs, D, YR,D, YRD, from 
which it tries to infer the pair of symbols bs, and bs, transmitted by Sı and S5, respectively. 
The derivation of the detector used by the destination, D, can be found in Section 4, and its 
performance analysis in Section 5. 


3. UEP—based network code design 


In (4), we have implicitly described the network code used by relay Ry with bg, = fr,(-,-). In 

this book chapter, four network codes (or NC scenarios) are investigated: 

e Scenario 1: bg, = bs, R, and bg, = bs. R- This network code corresponds to the working 
scenario in which relays R4 and Rz only decode-and-forward the signals received from 
sources 54 and 55, respectively (relay-only scenario) (Scaglione et al., 2006). 
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e Scenario 2: bg, = bs, R € ÎS R, and bg, = bs, R ® bs, Rə, Where © denotes bit-wise XOR 
operation. This network code corresponds to the working scenario in which relays R and 
Ry demodulate-network-code-and-forward the signals received from sources S; and Sp, 
respectively. Furthermore, they use conventional binary NC (XOR-only scenario) (Katti et 
al., 20082). 


e Scenario 3: bg, = bs, R ® bs, R, and bg, = bs. Ry This scenario corresponds to using a 
distributed network code obtained froma (4,2,2) UEP code (Van Gils, 1983, Table I), where 
a higher diversity gain has to be assigned to source 55. 


e Scenario 4: bg, = bs, R, and bg, = bs, R ® bs, R,. This scenario corresponds to using a 
distributed network code obtained froma (4,2,2) UEP code (Van Gils, 1983, Table I), where 
a higher diversity gain has to be assigned to source $1. 
Scenario 1 and Scenario 2 correspond to state-of-the-art distributed coding techniques 
(Rebelatto et al., 2010b), while Scenario 3 and Scenario 4 are the flexible network codes we 
are interested in studying in this book chapter. The reason why UEP coding theory can be a 
suitable tool to design distributed network codes for application scenarios in which different 
sources require a different diversity gain (see also Section 1) has its information-theoretic 
foundation in (Zhang, 2008). In fact, in (Zhang, 2008) it is shown that the minimum distance 
of a network code plays the same role as it plays in classical coding theory. Furthermore, 
from classical coding theory we know that the minimum distance of a linear block code 
directly determines the diversity gain over fully-interleaved fading channels (Proakis, 2000, 
Ch. 8), (Simon & Alouini, 2000, Ch. 12). In UEP linear codes, each systematic bit has its own 
minimum distance, and the set of these distances is known as separation vector (Masnick & 
Wolf, 1967), (Boyarinov & Katsman, 1981). From (Van Gils, 1983, Table I), the network code 
of Scenario 3 is a UEP distributed code with separation vector [2,3], which means that the 
minimum distance for the bits sent by S4 and Sz is equal to 2 and 3, respectively. Likewise, the 
network code of Scenario 4 is a UEP distributed code with separation vector [3, 2], which means 
that the minimum distance for the bits sent by S; and Sp is equal to 3 and 2, respectively. Thus, 
from (Zhang, 2008) it follows that by using a network code constructed from UEP coding 
theory we can individually assign different diversity gains to different sources. Also, note 
that, unlike (Xiao & Skoglund, 2009a), (Rebelatto et al., 2010b), this is obtained by neither 
using a non-binary Galois field nor introducing extra delays. The complexity and decoding 
latency of all the network codes studied in this book chapter are, on the other hand, the 
same. The downside is that only one source can achieve full-diversity. To the best of the 
authors knowledge, the adoption of UEP coding theory to design distributed network codes 
for multi-hop/cooperative networks with noisy and faded source-to-relay channels has never 
been addressed in the literature. Finally, we note that although, for analytical tractability and 
space constraints, only the two-source two-relay network topology is here investigated, UEP 
coding theory can be applied to generic multi-source multi-relay networks, by using, e.g., the 
codes available in (Van Gils, 1983, Table I). 


4. Receiver design 


We consider that the destination, D, uses a detector based on a low-complexity 
implementation of the Maximum Likelihood Sequence Estimation (MLSE) criterion. In 
particular, according to the MLSE criterion, given ys, p, Ys,p, YR,D, and yg,p introduced in 
Section 2, the proposed receiver estimates the distributed codeword that has most probably 


www.intechopen.com 


Flexible Network Codes Design for Cooperative Diversity 125 


been transmitted (Proakis, 2000). However, in order to keep the complexity of the detector 
at a low level, we introduce some simplifications in the analytical development with the aim 
of reducing the computational complexity and the a priori CSI required at the destination 
for optimal decoding. Of course, this leads to a sub-optimal receiver design, but allows the 
destination not to estimate the wireless channel over all the wireless links of the network. 

Before proceeding with the development of the detector, we need to identify the codebook, 
i.e., the set of distributed codewords that can be received by the destination in a noise-less and 
fading-less scenario. The codebook is denoted by C — f (1), c2), 2), 4) }, where c) is the 
(7) 


i 


is the i-th element of c) for i = 1,2,3,4. More specifically, in a 
G) G) (j) (j) 


noise-less and fading-less scenario we have: cï’ = bg, C5 = bs,, C4. = br and cù = bg,. 
As far as the NC scenarios described in Section 3 are concerned, the following codebooks can 
be obtained: 


e C = {0000, 0101, 1010, 1111} for Scenario 1 
e C = {0000, 0111, 1011, 1100} for Scenario 2 
e C = {0000, 0111, 1010, 1101) for Scenario 3 
e C = {0000, 0101, 1011, 1110} for Scenario 4 


j-th codeword of C, and c 


Detection at the destination encompasses two main steps, which involve physical and 
network layers, respectively. 


1. Step 1 (Physical Layer): From the received signals, ys, p, Vs, p, Vn, p, and yg, p, hard-decision 
estimates of [bs,,bs,,br,,br,| are provided by using a ML-optimum receiver, which 
exploits channel information on the source-to-destination and relay-to-destination wireless 
links (t = 1,2 and g = 1,2): 


bsp = argmin{|ys,p — VEmhs,p(1 — 2bs)?] 


bs,€{0,1} (5) 
bR,D = argmin{|yR,D — V Enhg, p(1 = 2bg )P] 
br,e{0,1} 
2. Step 2 (Network Layer): The hard-decision estimates ê = [(1,00,¢3,C4] = 


[bs, D: bs, p,b RiD,Ug,p] are input to the network layer, which uses a MLSE-optimum 
criterion to jointly estimate the bits emitted by the sources 5, and S», as follows (Proakis, 
2000), (Simon & Alouini, 2000): 


[5073 = arg max {Pr { 9s, bs | bs s Psp Pn i ôr Y} (6) 
bs, € (0,1), bs, c (0,1) 


where Pr {-} denotes probability. 


Since the network codes studied in Section 3 can be regarded as systematic linear block 
codes (Proakis, 2000), (6) can be re-written as follows: 


is, bo] = [P.e] args. {eff 


i E) \ ec 


Cs, Cy; 3 } (7) 
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Finally, by applying the Bayes theorem, by taking into account that the symbols emitted 
by the sources (and, thus, the codewords) are equiprobable, and by noticing that the 
hard-decisions € = [¢,€2,¢3,¢4] = |bs, p, bs, p, bn, p, Un,p] are independently computed, 


(7) simplifies as follows: 
Py 


[Bs] = [eP] = argmax mn I fe 
where, according to Section 3, we have: cl") = fh n J J ) and ol) = fn (a ? 2) ) : 


(PaPe Le 


By carefully looking at (8), we notice that the computation of Pr (à C (i) } fori = 1,2,3,4 


and j = 1,2,3,4 requires the knowledge of the channel gains over all the wireless links of 
the network, including the source-to-relay links to which the destination has not direct 
access. The availability of this information typically requires some cross-layer interactions 
between physical and network layers, along with the design of a so-called channel-aware 
detector (Chamberland & Veeravalli, 2007), (Di Renzo et al., 2009). With the aim to 
simplify the complexity of the detector, we retain two main assumptions (A) in this book 
chapter: A1) we consider that the destination has no knowledge of the channels over the 
source-to-relay links, and A2) we consider that the network layer only knows the fading 
distribution of the source-to-destination and relay-to-destination wireless links, but it 
does not know the exact realization of the channel gains. 


From A1), it can be readily proved that Pr fe 


ol) } fori = 1,2,3,4 and j = 1,2,3,4 
follows a Bernoulli distribution (Proakis, 2000) that does not depend on the source-to-relay 
links, as shown below: 


PP" 
Pr 4 ĉ o = (1- Psp)" a) pil) 


rele P | = a rao) tO n 


: (9) 
Pr {esc =a- Pg,p) ^U) PRU 


! "NS 
Pr E (LL Pg 5) A7) pall) 


where: i) Pxp = Q (v xo) with X = (5,, S2, R1, R2} is the error probability over 
the link from node X to node D, ii) 7 = 2E; / No, iii) Q(x) = (1/27) i exp( —? /2)dt 


is the Gaussian Q-function, and iv) d; (j) = Ap (ae) = le — V) for i = 1,2,3,4, 


1 


where dy (a,b) denotes the Hamming distance between two bits a and b (Proakis, 2000). 


From A2), as the network layer has no access to the actual fading gains but only knows 
their distribution, we need to replace Pxp in (9) with its average value, as follows: 
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Pxp = E(Pxp) = [ o (VIE) y, p E) 48 (10) 


where E{-} denotes the expectation operator computed over fading channel statistics, 
and Bino (-) is the Probability Density Function (PDF) of the fading power gain |hyp|*. 


Since, according to Section 2, we consider ii.d. Rayleigh fading channels, Pxp can 
be readily computed in closed-form as follows (Simon & Alouini, 2000) (for X — 


(51, S2, Ry, R2): 
AZ 1 202* 
pol. cT. I—9- 11 
merum V c 4 


Finally, by substituting (9) and (11) into (8), and computing the logarithm, we obtain, after 
some algebra, the result as follows: 


RA = Pur —  argmax fin (=) Tu a} 


i4 lee =? 
12 
$ a o (12) 
= argmax j In( 755) L 6-6 
fea lee px 


By taking into account that In (P/(1 — P)) is a negative (for P € [0,1/2]) factor that does 
not effect the outcome of the detector, (12) simplifies as follows: 


RA = Pur —  argmin b 4-4 (13) 
1,45 bec ES 


which turns out to be a (distributed) Minimum Distance Decoder (MDD) receiver for 
network decoding. 


In (13), it is important to note that multiple codewords of a codebook might have the same 
Hamming distance from the hard-decisions € = |64,62,65,64] = lbs, D: bs, Di br, D: br D] 
provided by the physical layer. In this case, we simply assume that the detector randomly 
chooses one of them with equal probability. 


5. Performance analysis 


The objective of this section is to develop an accurate analytical framework to compute the 
ABEP of the MDD receiver in (13) for the four NC scenarios described in Section 3. To this end, 
we compute first the BEP conditioned upon fading channel statistics, and then the average 
over the wireless channel. 


5.1 Conditional bit error probability (BEP) 
To compute the BEP, we take into account that: 
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1. For analytical tractability, we use union bound methods that require the estimation of the 
Pairwise Error Probability (PEP) for each pair of codewords of the codebook (Proakis, 
2000). 


2. We assume that the codewords of the distributed network code are equiprobable. 


3. We separately compute the BEP of sources 5, and S», since, as mentioned in Section 1 and 
Section 3, we are interested in showing that UEP-based distributed network codes provide 
different performance for different sources, according to the specified separation vector. 


Accordingly, the BEP of source S; for t = 1,2 can be upper-bounded as follows (t = 1,2)!: 


M ME . Ll 
BEP(5) = i Y Y PEP(S#) (h> jz) (14) 
h-ljf£hn-1l 


where PEP (5) (ji > j2) is defined as: 


PEP(5) (ji > jo) = Pr [eo — c2 


cH) T cen \ (15) 


where Pr} ct) — cU» lee Æ m = P { eli) E co») if cU) F cn and 


Pr { c(h) — «e» cu) x p = 0 if cU) = cun and it allows us to compute the BEP 
of each source individually. As a matter of fact, the detector might be wrong in estimating 
the transmitted codeword, but this does not necessarily lead to a decoding error in the bits 
transmitted by both sources. As an example, let us consider Scenario 3 and the transmission of 
c2) = 0111. If the receiver (wrongly) decodes e?) = c(4 = 1101, then this results in an error 
only for the bit emitted by 54, while there is no error for S». 


From (15) and by grouping together common PEPs, it can be shown that the BEP of 54 and 55 
in (14) simplifies as follows: 
BEP(5U = PEP(S) (1 = 3) + PEP($U (1 — 4) 
S S S (16) 
BEP(S2) = PEP(S (1 2 2) + PEP(%) (1 — 4) 


where we conclude that only three PEPs need to be computed to estimate the error 
performance, as it can be shown that: PEP(1— 2) = PEP‘) (1 > 2) = PEP 2) (1 > 2), 
PEP (1 => 3) = PEP‘) (1 > 3) = PEP‘) (1 3), and PEP (1 > 4) = PEP“) (1 4) = 
PEP(S2) (1 > 4). 


1 Note that, to simplify the notation, we avoid to emphasize that BEP and PEP are conditioned upon the 
fading channel. Instead, we use ABEP and APEP to denote the same functions when averaged over 
fading channel statistics. 
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The PEPs in (16) can be computed from (13). In particular, by direct inspection of (13), the 
generic PEP can be explicitly written as follows: 


cH) 4 cn \ 


EE { Di, > Dj, ie # cn \ F PP {Di, a Dj, [o E cn \ 


PEP(*) (ji > jo) =Pr (eo — c» 


(17) 


4 
where we have defined D; = 7 4 


ĉi — cÜ | for j — 1,2,3,4. Note that the second addend 


in the second line of (17) is due to the closing comment made in Section 4, where we have 
remarked that the detector randomly chooses with equal probability (i.e., 1/2) one of the two 
decision metrics D;, and Dj, in (17) if they are exactly the same. 


Let us now introduce the random variable: 


4 
Dj, - Dj, = 1 | 


D 


| (18) 


juj T 
Then, by denoting the probability density function of D;, ;, conditioned upon ci y Æ cl) 2) in 
(18) by gp, ;. (. cf 1) A ci 2) a the PEP in (17) can be formally re-written as follows: 


, : too ; ; 1 por ; 
PER) o ja) =f” soy, (e| fate f sona (e| He) ae a» 
Closed-form expressions of PEP(*) (ji — j2) are computed in Section 5.3. 


5.2 Average bit error probability (ABEP) 
The ABEP can be readily computed from (14) by exploiting the linearity property of the 
expectation operator. In formulas, we have: 


4 4 
ABEPÓ) = E { BEP(S)) } = ly: ys APEPÓ)(j >j) (20) 
4 i-i aei 


where APEP(S) (jj > j2) = E { PEP) Uh > jo). 


The APEPs in (20) can be computed by taking the expectation of (19) after computing the 
integrals. Closed-form expressions of these APEPs are given in Section 5.3. 


5.3 Average pairwise error probability (APEP) 
The closed-form computation of the APEPs in (20) requires the knowledge of the probability 


J 
opposed to many state-of-the-art research works, our system setup accounts for errors over the 


source-to-relay links. More specifically, (3) shows that the relays might incorrectly demodulate 
the bits transmitted by the sources. Even though the MDD receiver in (13) is unaware of these 
decoding errors, as explained in Section 4, they affect its performance and need to be carefully 
taken into account for computing the APEPs. 


density function gp, ,, ( | z 8T) in (19). In Section 2, we have mentioned that, as 
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More specifically, in Section 4 we have shown that the relays operate ina D-NC-F mode, which 
means that they perform two error-prone operations: i) they use the DF protocol for relaying 
the received symbols, and ii) they combine the symbols received from the sources by using 
NC. The accurate computation of the APEPs in (20) requires that the error propagation caused 
by DF and NC operations at the relays are accurately quantified. 


5.3.1 DF and NC operations: The effect of realistic source-to-relay channels 

As far as DF is concerned, the error propagation of this relay protocol in two-hop relay 
networks has already been quantified in the literature. In particular, in (Hasna & Alouini, 
2003) the following result is available. 

Given a two-hop, source-to-relay-to-destination (S-R-D), wireless network, the end-to-end 
(i.e., at destination D) probability of error, Psrp, is given by: 


Psrp = Psr + Prp — 2PsgPnp (21) 


where Psg and Prp are the error probabilities over the source-to-relay and relay-to-destination 
links, respectively. 
By taking into account the analysis in Section 4, it can be readily proved that Pop = 


Q ( Y ["3 and Pgp = Q (v d iol). The average end-to-end probability of error, 
Psrp, can be computed from (10) and (11), and by taking into account that channel fading 
over the two links is uncorrelated. The final result from (21) is: 


Psrp = E{Psrp} = Psr + Prp — 2PsrPRp = 2P — 2P? (22) 


Let us now consider the error propagation effect due to NC operations and caused by errors 
over the source-to-relay channels. In this book chapter, NC, when performed by the relays, 
only foresees binary XOR operations (see Section 3). Thus, we analyze the error propagation 
effect in this case only. The result is summarized in Proposition 1. 


Proposition 1. Let bs, and bg, be the bits emitted by two sources Sı and S» (see, e.g., (1)). 
Furthermore, let bs, and bs, be the bits estimated at relay R (see, e.g., (3)) after propagation through 
the wireless links S1-to-R and S5-to-R, respectively. Finally, let bg = bs, © bs, be the network-coded 
bit computed by the relay R. Then, the probability, Pg, that the network-coded bit, bg, is wrong due to 
fading and noise over the source-to-relay channels is as follows: 


Pg — Pr { (és, ® bs,) 7 (bs, ® bs,)} = Psp + Psr — 2PsirPs rp (23) 
where Ps, g and Ps, p are the error probabilities over the S4-to- R and S-to-R wireless links, respectively. 


Similar to the analysis of the DF relay protocol, it can be readily proved that Psr = Q (v "y lis) 
= 2 
and Par = Q (ys ) 


Proof: The result in (23) can be proved by analyzing all the error events related to the 
estimation of bs, and bs, at relay R. In particular, four events have to be analyzed: (a) no 
decoding errors over the S;-to-R and 5»-to-R links, i.e., bs, — bs, and bs, = bs,; (b) decoding 
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(a) No decoding errors (b) Decoding errors over the S; — R link 


bs, | bs, | bs, @bs, | bs, | bs, | br | — | bs, | bs, | bs, Sbs, || bs, | bs, | br | 
0 


S2 S2 S2 

0 0 0 0 0 1 
1 0 1 1 

0 1 0 1 

1 1 1 0 


(d) Decoding errors over both links 


bs, | bs, | bs, Obs, | bs, | bs, | br | 


e e 


| bs, | bs, | bs, Obs, || bs, | 
0 0 0 


Sı 
0 


0 0 0 1 1 0 
1 
1 
0 


Table 1. Error propagation effect due to NC at the relays for realistic source-to-relay channels. 


errors only over the 5,-to-R link, i.e., bs, # bs, and bs, = bs,; (c) decoding errors only over the 
S5-to-R link, i.e., b sı = bs, and bs, # bs,; and (d) decoding error over both 5,-to-R and S-to-R 
links, i.e., bs, # bs, and bs, # bs,. These events are summarized in Table 1. In particular, we 
notice that errors occur if and only if there is a decoding error over a single wireless link. 
On the other hand, if errors occur in both links they cancel out and there is no error in the 
network-coded bit. Accordingly, Pr can be formally written as follows: 


Pg = Pr { bs, Æ bs, } + Pr (bs, Æ bs, } —2Pr ls, Æ bs, and bs, Æ bs, } 


— Pr fbs, Æ bs, } “pure fbs, Æ bs, } =2Pr fbs, £ bs, } Pr fbs, Æ bs, } d 


which leads to the final result in (23). This concludes the proof of Proposition 1. 


Finally, we note that, from (23), the average probability of error at the relay with NC, Pr, can 
be computed from (10) and (11), and by taking into account that the fading over two links is 
uncorrelated. The final result from (23) is: 


Pg = E {PR} = Pag + Ps,g — 2Ps RPs a = 2P — 2P? (25) 


Very interestingly, by comparing (22) and (25) we notice that DF and NC produce the same 
error propagation effect. Thus, by combining them, as the network codes in Section 3 foresee, 
we can expect an error concatenation problem. In particular, by combining the results in (22) 
and (25), the end-to-end error probability of the bits emitted by sources S4 and $5 and received 
by destination D (denoted by Ps,(r,R,)p and Ps,(r,R,)p, respectively) can be computed as 
shown in (26)-(29) for Scenario 1, Scenario 2, Scenario 3, and Scenario 4, respectively: 


Rud = Ps rR, + PRD — 2Ps, R PRID (26) 


Ps, (RiR) D = PsiR + Prap — 2Ps gm; Pp 
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Ps (RiR2)D = [Ps,n, + Ps,r, — 2Ps,n, Ps;n;] + Priv — 2[Ps;n, + Ps,n, = 2Ps n, Ps,n;]Pnip 
Ps; R,R2)D — [Ps R; sr Pim a 2Ps, R, Poss] "E PRD = 2[Ps. R; + Ps, R, = 2Ps n Pss; | Pop 
(27) 
Psi RQR2)D — [Ps F Ps R; mi 2Ps, R Pss] T Priv a Z| Por: + Ps R, m 2Ps, R Ps, a, | PRD 
Ps, (RR,jD = Psn, + Prod — 2Ps,R,PR»D 
(28) 


Fl RiR2)D = Ps,n, + Priv — 2Psg, Pup 

Ps,(RiR2)D = Bs + Ps,n, — 2Ps, n, Ps,n;] + Pn;p — 2[Ps n, + Ps,n, — 2Ps, n, Ps, n, |Pn;p 
Q9) 

The average values of Ps (ggp and Ps (RiR)D Les Ps(ngp = E [Psp] and 

Ps Reo = E LPs. RiR2) p] can be computed by using arguments similar to (22) and (25). 

The final result is here omitted due to space constraints and to avoid redundancy. 


5.3.2 Closed—form expressions of APEPs 

From (16) and (20), it follows that only three APEPs need to be computed, for each NC scenario 
in Section 3, to estimate the ABEP of both sources. Due to space constraints, we avoid to 
report the details of the derivation of each APEP for all the NC scenarios. However, since 
the derivations are very similar, we summarize in Appendix A the detailed computation of a 
generic APEP. All the other APEPs can be derived by following the same procedure. 

In particular, by using the development in Appendix A the following results can be obtained: 
Scenario 1: 


APEP(1—2) = PspPs(mnjp-(1/2)(1— Psp) Ps,mnjp 
+ (1/2) (1 - Ps. np) Pen 


APEP(1—3) = Ps DPs (RR) + (1/2) (1 — Psp) Psi (RiR)D 
+ (1/2) (1 E Ps, (RiRa)D) Psp 


APEP(1—4) = (1/2) (1— Ps,p) (1 — Ps,p) Ps, (mn; ypPs,(R;R;)D 
(1/2) (1.— Psp) (1 — Ps, (n g;jp ) Ps;pPs, (n Rj)p 
(1/2) (1 — Psp) (1 — Psi np) PspPs (nn) (30) 
(1/2) (1 = Psp) (1 — Ps (RiR)D ) Ps;pPs,(n;)D 
(1/2) (1— Psp) (1 — Psi np ) Ps;pPs (gg) 
(1/2) (1 — Ps, (RRID p) (1 " Ps. np) Pon Pep 


de » D) Ps,pPs, (m n;)DPs, (m Rj)D 

(1 — Ps, p) Ps; DPs (gn; jpPs,(R R)D 
Le Py Rikin ) Ps,pPs,pPs, (mn; )p 
1 — Psi n;jp ) PspPs,pPs, (n gp 
Ps DPs, DPs, (n g;)DPs, (Ri R;)D 


++ +++ é k 4 4 4+ 4 
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Scenario 2: 

APEP(1— 2) = (1-Psp)Ps(nnj)pPs,(n RD 
ct. (T= Ps (RRD ) PspPs,(R RD 
+ (1-Psmnjp) PspPs(ngj)p 
+  PspPs(nnpPs,(nmj)D 

APEP(1— 3) = (1-Psp)Ps(nnjpPs,(nnj)D (31) 
an He P, S, (R4 R2)D ) Ps,pPs,(R,Ro)D 
+ (1-Psqmnyp) PspPs(n gp 
+ Ps DPs (n gjpPs,(n RD 

APEP(1—4) = (1/2)(1— Psp) Psp + (1/2) (1 — Psp) Psp 
+ PspPs,D 

Scenario 3: 

APEP (1> 2) = (1- Psp)Ps(nnjpPs(nn)p 
ae qe P, S. (Ri R85)D ) Ps,pPs, (s RD 
+ [1-Psgmnjp)PspPs(nR)p 
+ PspPs n g;ypPs, (gj Rj)D 

APEP(1—3) =  PspPs(mmjp- (1/2) (1 — Ps,p) Psi(RiR)D 

: : (32) 

+ (1/2) (1 " Ps. RyRa)D) Psp 

APEP(1—4) = (1-Psp)Ps(nnjpPs(nR)p 
+. (1- Ps (RRD ) PsiDPs,(RiR)D 
+ = Ps,(R,Ro)D ) laD RRD 
23 bon. Ri R;)DPS (Ri R5)D 
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Scenario 4: 

APEP(1—52) = (1/2) (1—Ps,p) Ps,mn;p + Ps,pPs,(mn;)p 
+ (1/2) (1 — DT) Psp 

APEP(1—3) = (1- Psp) Ps,(R,R»)DPs)(R,Ro)D 
+ dl-Psmnyp)PspPsnRjp 
+ [i-Psi(injp)PspPs n np (33) 
+ PspPs,(nn;pPs,(R; R;)D 

APEP (1 — 4) = ( — Psp ) Ps,pPs (n R,)D 
+ (L—Ps5)PspPstgRD 
+ (1- Paran) PspPs,p 
+ Pe prop seein 


5.4 Diversity analysis 

Let us now study the performance (ABEP œ) of the MDD receiver for high SNRs, which allows 
us to understand the diversity gain provided by the network codes described in Section 3 
(Wang & Giannakis, 2003). To this end, we need to first provide a closed-form expression of 
the ABEP of Sı and Sz from the APEPs computed in Section 5.3.2. By taking into account 
that the wireless links are i.i.d. and that the average error probability over a single-hop link is 
given by P in (11), from (20), (30)-(33), and some algebra, the ABEPs for Scenario 1, Scenario 2, 
Scenario 3, and Scenario 4 are as follows, respectively: 


Scenario 1: ABEP) = ABEP») =  (1/2)P, + (1/2)P3 + (1/2)P, P5 + (1/2)P, P + 
(1/2),P, + (1/2BbP, + (1/2bPB, — (1/2b,bP, — (1/2)P,P),P, — 
(1/2)P4 BP, — (1/2) SP. — (1/2)P1 P5 P4 P4, where we have defined P; = P) = P 
and P = P4, = 2P — 2P?. 

Scenario 2: ABEP(®) = ABEP?) = (1/2)P,; + (1/2)P + P,P; + P,P, + DP, — PPP, — 
P, P3P,, where we have defined P, = P; = P and P4 = P, = 3P — 6P? + 4P°. 

Scenario 3: ABEP) = (1/2)P, + (1/2)P3 + PP; + P,P, + P,P, — P PaP, and ABEP(2 = 
P,P) + P Ps + P4P4 + 2P5 P4 + P3P4 — 2P, PP, — 2P>P3P,, where we have defined 
P, = P, = P, P = 3P — 6P? + AP?, and Py = 2P — 2P?. 

Scenario 4: ABEP) = P,P, + 2P,P; + P,P, + LB + P3P, — 2P,P5 P, — 2P PP, and 
ABEP(?) = (1/2)P + (1/2)Py + P,P) + P,P; + 2P5 P4 — PaP, — 2P4P5 P5, where 
we have defined P, = P; = P, P = 2P — 2P?, and P4 = 3P — 6P* + 4P°. 


From the results above, we notice that in Scenario 1 and Scenario 2 both sources have the 
same ABEP. Furthermore, for all the NC scenarios we can easily compute ABEP% and the 
diversity gain (Div) of 54 and S5, as shown in Table 2. In particular, from Table 2 we observe 
that, by using UEP coding theory for network code design (i.e., Scenario 3 and Scenario 4), at 
least one source can achieve a diversity gain greater than that obtained by using relay-only 
or XOR-only network codes (i.e., Scenario 1 and Scenario 2). Furthermore, this performance 
improvement is obtained by increasing neither the Galois field nor the number of time-slots 
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Scenario 1 | (3/2)P (3/2)P 
P 


Scenario 2 P 
Scenario 3 2P 16P? 
Scenario 4 16P? 2P 


Scenario 1 


Scenario 2 


Scenario 3 


Scenario 4 


Table 3. ABEP% of S; and Sz and diversity gain with ideal source-to-relay channels. 


(Rebelatto et al., 2010b). Finally, by studying the diversity gain provided by the network 
codes obtained from UEP coding theory in terms of separation vector (SP), we observe that 
the achievable diversity gain is equal to Div — SP — 1. From the theory of linear block codes, 
we know that this is the best achievable diversity for a (4,2,2) UEP-based code that uses a MDD 
receiver design at the destination (Proakis, 2000), (Simon & Alouini, 2000). Better performance 
can only be achieved by using a more complicated receiver design, which, e.g., exploits CSI at 
the network layer. 


5.5 Effect of realistic source-to-relay channels 

In Section 2, we have mentioned that the relays simply D-NC-F the received bits even though 
the source-to-relay channels are error-prone, and so the transmission is affected by the error 
propagation problem. Thus, it is worth being analyzed whether this error propagation effect 
can decrease the diversity gain achieved by the MDD receiver or whether only a worse coding 
gain can be expected. To understand this issue, in this section we study the performance of 
an idealized working scenario in which it is assumed that there are no decoding errors at the 
relays. In other words, we assume bs, RIZ bs, for t = 1,2 and q = 1,2 in (3). In this case, 
the expression of the ABEP for high SNRs can still be computed from (20) and (30)-(33), but 
by taking into account that P = Psp = Ps,p = Ps (RiR) D = Ps,(R,R»)D- The final result of 
ABEP for S4 and S5 is summarized in Table 3. 

By carefully comparing Table 2 and Table 3, we notice that there is no loss in the diversity gain 
due to decoding errors at the relay. However, for realistic source-to-relay channels the ABEP 
is, in general, slightly worse. Interestingly, we notice that Scenario 2 is the most robust to error 
propagation, and, asymptotically, there is no performance degradation. 


6. Numerical examples 


In this section, we show some numerical results to substantiate claims and analytical 
derivations. A detailed description of the simulation setup can be found in Section 2. In 
particular, we assume: i) BPSK modulation, ii) o = 1, and iii) according to Section 5.5, both 
scenarios with and without errors on the source-to-relay wireless links are studied. 
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Scenario 1 Scenario 2 
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Fig. 2. ABEP against E; / No. Solid lines show the analytical model and markers Monte Carlo 
simulations (0% = 1). 


The results are shown in Figure 2 and Figure 3 for realistic source-to-relay links, and in 
Figure 4 and Figure 5 for ideal source-to-relay links, respectively. By carefully analyzing 
these numerical examples, the following conclusions can be drawn: i) our analytical model 
overlaps with Monte Carlo simulations, thus confirming our findings in terms of achievable 
performance and diversity analysis; ii) as expected, it can be noticed that the ABEP gets 
slighlty worse in the presence of errors on the source-to-relay wireless links for Scenario 1, 
Scenario 3, and Scenario 4, while, as predicted in Table 3, the XOR-only network code (Scenario 
2) is very robust to error propagation and there is no performance difference between Figure 
2 and Figure 4; and iii) the network code design based on UEP coding theory allows the MDD 
receiver to achieve, for at least one source, a higher diversity gain than conventional relaying 
and NC methods, and without the need to use either additional time-slots or non-binary 
operations. 

More specifically, the complexity of UEP-based network code design is the same as relay-only 
and XOR-only cooperative methods. For example, by looking at the results in Figure 3 and 
Figure 5, we observe that the network code in Scenario 3 is the best choice when the data sent 
by Sz needs to be delivered i) either with the same transmit power but with better QoS or ii) 
with the same QoS but with less transmit power if compared to S4. The working principle 
of the network code in Scenario 3 has a simple interpretation: if S2 is the "golden user”, then 
we should dedicate one relay to only forward its data without performing NC on the data of 
$1. A similar comment can be made about Scenario 4 if S4 is the "golden user". This result 
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Fig. 3. ABEP against Em / No. Solid lines show the analytical model and markers Monte Carlo 
simulations (o5 = 1). 


highlights that, from the network optimization point of view, there might be an optimal choice 
of the relay nodes that should perform relay-only and NC coding operations. By constraining 
the relays to perform simple operations (e.g., to work in a binary Galois field), this hybrid 
solution might provide better performance than scenarios where all the nodes perform NC. 
However, analysis and numerical results shown in this book chapter have also highlighted 
some important limitations of the MDD receiver. As a matter of fact, with conventional 
relaying and NC methods only diversity equal to one can be obtained, while with UEP-based 
NC at least one user can achieve diversity gain equal to two. However, the network topology 
studied in Figure 1 would allow each source to achieve a diversity gain equal to three, as 
three copies of the messages sent by both sources are available at the destination after four 
time-slots. This limitation is mainly due to the adopted detector, which does not exploit 
channel knowledge at the network layer and does not account for the error propagation 
caused by realistic source-to-relay wireless links. The development of more advanced 
channel-aware receiver designs is our ongoing research activity. 


7. Conclusion 


In this book chapter, we have proposed UEP coding theory for the flexible design of network 
codes for multi-source multi-relay cooperative networks. The main advantage of the proposed 
method with respect to state-of-the-art solutions is the possibility of assigning the diversity 
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Fig. 4. ABEP against Em / No. Solid lines show the analytical model and markers Monte Carlo 
simulations (o5 = 1). Ideal source-to-relay channels. 


gain of each user individually. This offers a great flexibility for the efficient design of 
network codes for cooperative networks, as energy consumption, performance, number of 
time-slots required to achieve the desired diversity gain, and complexity at the relay nodes 
for performing NC can be traded-off by taking into account the specific and actual needs of 
each source, and without the constraint of over-engineering (e.g., working in a larger Galois 
field or using more time-slots than actually required) the system according to the needs of the 
source requesting the highest diversity gain. 

Ongoing research is now concerned with the development of more robust receiver schemes at 
the destination, with the aim of better exploiting the diversity gain provided by the UEP-based 
network code design. 
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Fig. 5. ABEP against Em / No. Solid lines show the analytical model and markers Monte Carlo 
simulations (o5 = 1). Ideal source-to-relay channels. 
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A. Appendix: Proof of (30)-(33) 


To understand how (30)-(33) are computed, in this section we provide a step-by-step 
derivation of the computation of APEP(1— 3) for source S, and Scenario 1, i.e., 
APEP!) (1 + 3) = Pr {0000 — 1010}. Note that since ci Ü xx ce 1, we avoid to 
emphasize, for the sake of simplicity, this conditioning in what follows. Other APEPs, for all 
the other scenarios, can be obtained with a similar analytical derivation. 

From (19), the PEP can be computed as follows: 


+00 0* 
PEPO 5 3)- [eo (dte 5 |. sos (dt (04) 


where gp, , (-) is the probability density function of random variable D; 3: 


Dis - Y [à - «?|- |a - |] = (35) 


i=] i=1 


where guo = 


A 1 
6 —d |= 


6 — ej | fori = 1,2,3,4. 


By direct inspection, it is possible to show that guo» fori = 1,2,3,4 are independent Bernoulli 


distributed random variables with probability density function as follows: 


28, (£) = (L— Psp)? (E +1) + Papé(o — 1) 

8g, (E) = o (C " 
88, (6) = (1 = Psi Rp) ô (6 +1) + Ps, (agp? (6 — 1) PR 
85, (E) = 8 (E) 


where ô (-) is the Dirac delta function. 
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It is relevant to notice that gg, (¢) = gg, (6) = ô (Ẹ), ie., B2 = B4 = 0 with unit probability, 
a) _ (3) (1) _ (3) 


because c) = c} and c} = c; ,and, so, regardless of the estimates €? and ĉ4 provided by 
the physical layer, we always have |é2 — ci | — |e) — ce | = 0 and 6&4 — | = 64 = ce | =Q. 


Since geo» fori — 1,2,3, 4 areindependent random variables, the probability density function 
of D, 5 in (35) can be computed via the convolution operator: 


8n,4 (6) = (sp, B gp, G 95, & gg, ) (6) — (sp, ® gp, ) (¢) 
= ja — P p) Psi (RıR)D + (1 E. Psi (RiRa)D ) Psp| ô (£) (37) 
+(1= Psp) e = Psi (RiRa)D ] ô (6 +2) + PsiDPsi (RRD (6 — 2) 


where & denotes convolution. 
Furthermore, by substituting (37) into (34) we can get the final result for the PEP: 


PEP) (1 — 3) = (1/2) (1 — Psp) Ps no + (1/2) (1 — Ps, rks) Psi 


+ Ps. DPs, (RiR) D 


(38) 
Finally, the APEP can be computed by simply taking the expectation of (38) and by considering 
that fading over all the wireless links is independent distributed: 

APEP() (1 - 3) = E {PEP (1 + 3)} 


= (1/2) (1 — Psp) Ps, gp + (1/2) (1 = Ps,(RiRa)D) Psp 69) 
+ Ps, pPs (Ry RD 


We observe that (39) coincides with (30), and this concludes our proof. 
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